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Abstract. For Gaussian initial conditions the perturba- 
tion theory predicts a very specific hierarchy for the pro- 
jected matter p-point correlation functions. In the small 
angle approximation and assuming a power-law spectrum 
I derive the exact expressions of the coefficients s p relat- 
ing the averaged p-order angular correlation function, U p 



V) to the second one, u>„ 



0> 



These results are valid 



for any selection function, but for a top-hat angular filter 
only. These coefficients are found to be significantly higher 
than their 3D counterparts, S p = £, p /~^2 

For the coefficient S3 I discussed the accuracy of the 



[T 1 small angle approximation by computing, for particular 
I examples, its angular dependence with Monte-Carlo nu- 
merical integrations. It is found that the accuracy of the 
small angle approximation for 6 ps 1° slightly depends 

q>^ on the selection function. Using the selection function ex- 

00 pected for galaxy catalogues the approximation is found 
to be reasonably good. 

The measurements of the s„ parameters made in the 
APM angular survey are found to give systematic lower 

q>^ values than the theoretical predictions. How significant 
this discrepancy is and what the implications would be for 
galaxy formation models is discussed in the last section. 
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1. Introduction 

The study of the statistical properties of the galaxy dis- 
tribution had been done initially in angular surveys such 
as the Lick catalogue (Shane & Wirtanen 1967, Seldner 
et al. 1977). The more recent apparition of redshift sur- 
veys led people to study in more details of the galaxy 3D 
distribution, and most of the results obtained in 2D cata- 
logues have been confirmed. In particular the index of the 



2-point correlation function initially determined by Lim- 
ber (1954), Neyman, Scott & Shane (1956) and Groth & 
Peebles (1977) has been reobtained by the analysis of the 
redshift surveys. 

However, due to the huge number of objects contained 
in angular surveys, they still allow the most precise sta- 
tistical studies. In particular one can more reliably mea- 
sure high-order correlation functions in 2D catalogues 
than in 3D. The APM Galaxy Survey with more than 
1,300,000 objects is at present the largest catalogue avail- 
able and Gaztanaga (1994) was able to measure the av- 
eraged galaxy angular p-point correlation functions up to 
p = 9. Similar results are beyond the present status of the 
observations for redshift surveys (Bouchet et al. 1993). 
Moreover in 2D catalogues there are no redshifts effects 
that can mix the density field and the velocity field prop- 
erties, which makes comparisons with theoretical predic- 
tions more difficult. The main difficulty brought by the 
analysis of the angular surveys is that a given angular 
scale corresponds to various physical scales. Limber (1954) 
pioneered the expected relationship between the angular 
2-point correlation function and the real space correlation 
function. This relation, and its approximation for small 
angles (eg., Peebles 1980, §50), proved to be efficient and 
accurate. 

In the strongly nonlinear regime, our theoretical knowl- 
edge of the quantitative behaviour of the density field is 
still weak. A-body simulations have been widely used to 
improve our understanding of this complicated part of 
the dynamics. The size of the catalogues involved in 2D 
measurements, however, still challenges the largest simu- 
lations. On the other hand, rapid progress have been made 
recently in the mildly nonlinear regime where perturbation 
theory can be applied. This led to a better understanding 
of the early departure from the linear behaviour. In par- 
ticular it has been shown that in such a regime, the mean 
value, £ , of the p-point correlation function £ p (xi, . . . ,x p ) 

in a spherical cell is proportional to £^ , 
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(1) 



vvnere Op are coeiiicieins uepenuing on une snape 01 
the power spectrum only (Fry 1984, Goroff et al. 1986, 
Juszkiewicz, Bouchet & Colombi 1993, Bernardeau 1994a, 
b). This hierarchy is specific of Gaussian initial conditions. 
For non-Gaussian initial conditions similar hierarchy may 
or may not be obtained. It seems to depend explicitly on 
the specificities of the model. Early studies (Jaffe 1994) 
tend to demonstrate that the hierarchy (1) is a robust 
result and holds even for mildly non-Gaussian initial con- 
ditions. Numerical studies show that the analytical results 
obtained in the quasilinear regime are very accurate to de- 
scribe the correlation hierarchy up to a = 1 (Juszkiewicz 
et al. 1994, Bernardeau 1994, Lokas et al. 1994, Baugh, 
Gaztanaga & Efstathiou 1994). The 3D statistics of the 
one-point density PDF in the quasilinear regime has thus 
been investigated in details, but for the reasons men- 
tionned before comparisons with the observations remain 
questionable. It may be worth to note that numerical sim- 
ulations show that the nonlinear corrections tend to in- 
crease the value of the coefficients S p (Bouchet & Hern- 
quist 1992, Colombi et al. 1994, Lucchin et al. 1994). 

However, so far, substantial progress has been made 
only for the statistical properties of the density field when 
it is averaged in a spherical cell. The properties of a 
smoothed angular density field corresponds rather to a 
filtering in a cone of a given solid angle. It implies that, 
for 2D catalogues, the connected part of the moments, 
the cumulants c p (for a definition, see Balian & Schaeffer 
1989), are expected to exhibit a similar specific behaviour 
(Schaeffer 1987), 

Cp — Sp C?j ; (2) 

but with different coefficients s p . These coefficient s p are 
expected to depend on the power spectrum as well as on 
the shape of the selection function, F(r). Note that the 
property (2) can be expressed in term of the p-point an- 
gular correlation function, w p (ji , . . . , j p ). The relation (2) 
implies that the average of the p-point angular correlation 
functions, U p , within a cell of angular radius 6o follow the 
hierarchy, 

LOp — Sp . (3) 

The hierarchy (3) is indeed found to describe accu- 
rately the statistics of both the galaxies and the galaxy 
clusters. The parameters s p are measurable up to p = 6 
for the clusters and up to p = 9 at small scales for the 
galaxies (Cappi & Maurogordato 1994, Gaztanaga 1994). 
This result has been used as an argument in favor of the 
gravitational instability scenario and for the absence of bi- 
ases in the galaxy distribution. Unfortunatly the parame- 
ters S p , corresponding to another smoothing scheme, can- 
not be simply derived from the observations in 2D cat- 
alogues. The approximation used by Toth, Hollosi and 
Szalay, (1989), and Gaztanaga, (1994), has been improp- 
erly used since it is valid only for exact tree models for 



une p-poim correlation luucuious. iiiereiore, iu snouiu nou 
be used for perturbation theory calculations. It gives how- 
ever an estimation of these coefficients that is attractively 
close to the measured values (Gaztanaga 1994, Frieman & 
Gaztanaga 1994). It is then worth to derive more exactly 
the quantitative behaviour of the high order cumulants 
from the perturbation theory. 

Throughout this paper, I make the hypothesis that the 
galaxies are good tracers of the matter field. In §2 I present 
the expected values of the coefficients s p as given by the 
perturbation theory in the small angle approximation. §3 
is devoted to a discussion on the accuracy of the small an- 
gle approximation on the value of S3 and the last section is 
devoted to a comparison of these results with observations 
and to comments. 

2. The projection effects in perturbation theory: 
the correlation hierarchy in the small angle 
approximation 

2.1. Physical assumptions 

The basic physical assumption is that the local number 
density of observable objects, n(x), is assumed to be pro- 
portional to the local matter density p(x) multiply by a 
selection function which is function of distance, F(r), 

n(x) = n F[\x\] p(x)/p , (4) 

where po is the mean density of the Universe. The quantity 
noF(r) is then the mean density of observable objects at 
the distance r. The selection function obviously depends 
on the sample that is considered. In the text I discuss the 
results for different shape of selection functions. 

The projected number density of objects, Co, per solid 
angle is 

/>oo 

a dtt = dtt n / r 2 drF(r), (5a) 
Jo 

and the local angular density, at the angular position de- 
fined by the unit vector 7, is 

/>oo 

a(f) = no / r 2 drF(r)p(r~f)/ po- (56) 
Jo 

The correlation properties of the angular density then 
depends on the hypothesis made for the matter density 
field. The complete series of the angular correlation func- 
tions, Up(7j, . . • ,7 p ), are difficult to measure. What can 
be more confidently measured are the spatial averages of 
these functions in a cell of angular size #0 , 

ZJ p = |d 2 7l W( 7l )...d 2 7p W( 1 p)ujp( ll) ... )lp ) (6) 

where W(f) is unity if the angle between 7 and a given 
direction 7 is less than #0 and zero otherwise (the val- 
ues of Up, according to the isotropic principle, should not 
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given by the cumulants of the number density probability 
distribution function when the angular field is filtered by 
a top-hat filter of radius 60 ■ 

The last important assumptions are that, at suffi- 
ciently large scale, the matter correlation properties can 
be described with the perturbation theory; and that the 
initial conditions were Gaussian. In most of the follow- 
ing calculations I will also assume that these scales are 
reached in the present 2D catalogues even for small fil- 
tering angles, so that for these calculations I will assume 
that 



< 1. 



(7) 



The validity of this assumption is discussed in §3. 

The assumption that perturbation theory can de- 
scribed the matter correlation properties implies that the 
three-point correlation function, for instance, has to be 
calculated using both the 1 st and the 2 nd order for the 
density field in perturbation theory. The properties of the 
linear density field are entirely given by the shape of the 
power spectrum P(k). It is defined in such a way that the 
Fourier transforms of the linear density contrast, given by 



bles 1980, §50) with such an approach. To do so I use the 
fact that, when ^Cl, 



J d 2 f W(f ) exp(i x k.f ) k, 
2-7T 9 2 exp[i x k cos^)] 



Ji[k x sin(flfc) 6> ] 
k x sin(#j;) #0 



(11) 



where k is the modulus offc, 0^ is the angle between k and 
7 and J\ is the Bessel function of the first kind. Using 
the relation (11) it is possible to derive the expression of 
the variance, U 2 , 



U> 2 



— — — — 2 / rJdnF(ri) / r 2 dr 2 F{r 2 ) 
[J r A drF(r)\ Jo Jo 



+00 



d k r f d 2 fcj_ 



, 2tt J (2tt) 2 
x W 2 (k L n0 o ) W 2 (k L r 2 6 ) 



P(k) exp[ifc r (ri - r 2 )] 



where I wrote 



k = k r f + fcj_ 



(12) 



(13) 



P (1) (x)/po 



d 3 fc 

(2tt) 3 / 2 



6(k) exTp(ikx), 



are Gaussian variables of moments 



^dirac (fei + k 2 )P(k 1 ). 



(8) 



In the following the power spectrum is approximated by 
a power law, 



P(k) = AV' 



(9) 



normalized to the present time. The index n can in princi- 
ple be determined from the shape of the angular two-point 
correlation function. 



2.2. The variance and the skewness 

To start with let me remind the expression of the density 

ad 



field at the 2 nd order 
f P\x) 



Po 



d 3 fcl d 3 fc 2 



(2tt) 3 / 2 (2tt) 3 / 
5 k\.k 2 
7 + ~fcF 



2 <5(fci) 6(k 2 ) exp[i(fci + k 2 )x] 



2 / fei.fea 

7 \hh 



(10) 



This expression is exact for an Einstein-de Sitter Universe 
only. This cosmological hypothesis is assumed in §2.2 and 
2.3. The dependence of the results with SI is discussed in 
§2.4. 

As the expression of the 2 nd order density field is sim- 
ple only in Fourier space we are led to derive the expression 
of the variance in this space. The first step is to rederive 



W 2 (x) = 2 



Ji(x) 



(14) 



In the limit of small angle #0 , the cell in which the integra- 
tion is made is extremely elongated so that the component 
of k along the line of sight is negligible compared to its 
orthogonal component, which implies 



k ki_. 



(15) 



With such an approximation the integrale (12) simplifies 
dramatically. Indeed the integration over r is straightfor- 
ward and implies that r\ = r 2 so that (12) reads, 

^= Ae i" s: S" iTF ' i :i r^-wiu). ,16) 



2tt 



[Jo°° r ' 2< irF(r)] Jo 



It is possible to show that the expression (16) is consistent 
with the Limber equation, 

/•OO /-oo 

w{9) (x r 4 F 2 (r)dr / £ 2 (V* 2 + r 2 6 2 ) dx, 



showing that (15) is indeed equivalent to the usual small 
angle approximation. 

The derivation of the skewness of the angular density 
field is based on similar approximations. It involves the 
integration over two different wave vectors, k\ and k 2 . 
Their components along 7 are then also negligible com- 
pared to their orthogonal components. As a result, in the 
expression of the 2 nd order density field the geometrical 
part 5/7 + ki .k 2 jk\ + 2/7 (fei .k 2 jk 2 k\) 2 involves only the 
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averaged three-point correlation function then reads, 



6.0. n iyn-ui uer cui i eiuiiuns 



_ 6A 2 6~ 2(n + 2) / °° r 8-2(n+3) drf 3 (r) 

x l\ +n W 2 (h) l\ +n W 2 (l 2 ) W 2 (\h + h\) 
x [b/7 + h.h/ll + 2/7 (h.l 2 /hh) 2 ] ■ 



d 2 h d 2 l 2 



(17) 



The calculation of the integrale over the vectors l\ and 
l 2 turns out to lead to simple integrations because of spe- 
cial geometrical properties of the function W 2 . Indeed we 
have 



W 2 (\h+l 2 \) [l-cos 2 (^)] = 7T W 2 (h) W 2 (l 2 )(l8) 



where <p is the angle between l\ and l 2 . We can also notice 
that 



2 dJ o (0 



/ dl 
and that 

/0 



-W 2 (l) 



,-2-k 

/ d^ J (|Zi + Z 2 |) = 2tt J (/i) J (/ 2 ) 
Jo 



which by differentiating with respect to l 2 leads to the 
relationship 



d^ W 2 (\h+l 2 \) [l + cos(p)/i// 2 ] 



= 2tt J (h) W 2 (l 2 ) 
= 2tt W 2 (l 2 ) W 2 (h) + 



(21) 



The properties (18) and (20) can be shown using similar 
decompositions (in sum of product of Bessel functions) 
than the ones used for the Fourier transform of the top-hat 
window function in 3D (Bernardeau 1994a). The demon- 
strations will not be given. 

Using the relationships (18) and (21) one can then cal- 
culate the expression (17). The result reads 



S3 



LU3 

= R3 



36 3 , 
T -2(" + 2 ) 



with 



/ °° r 8 - 2 ( n + 3 )drF 3 (r) / °° r 2 ArF(r) 
[/~r 5 -("+ 3 )drF 2 (r)] 2 



(22) 



(23) 



The expression of the measured coefficient S3 then re- 
duced to a simple expression where its dependence with 
the power spectrum index, and the shape of the selection 
function can be easily investigated. 



The generalization of the results obtained in part 2.2 to 
high-order correlation function is actually quite straight- 
forward. Let me remind the form of the general expression 
of the density contrast at the p th order (Goroff et al. 1986), 



P' 



.00 



Po 



d 3 fc! 

(2tt) 3 / 2 
x # p (fci, 



d 3 fc p 
(2tt) 3 / 2 
. , k p ) 



8(k 1 )...8(k p ) 



(24) 



where N p is a homogeneous function of the wave vectors 
k\ , . . . , kp . 

The p-order correlation function in the quasi-linear 
regime is given by product of the density field at various 
order following a kind of tree form (Fry 1984, Bernardeau 
1992). When one applies the approximation (15) for each 
of the wave vectors involved in the tree product the coef- 
ficient s p is expected to take the form (Toth et al. 1989, 
Gaztanaga 1994) 



(19) s p = R p T, p (n) 



where T, p is a function of n only and 



(20) R P 



M[~ 2 (3) M p [3p-(p-l)(n + 3)] 

MP-\3-n) 



with 



Mp(a) 



F p (r). 



(25) 



(26a) 



(266) 



This is a key property which implies that, in the limit of 
small angles, the dependence with the shape of the selec- 
tion function factorizes away and is very simple. Note that 
the parameters R p are generally greater but close to unity. 

The coefficients T, p are however not directly related 
to the S p coefficients. Similarly to those coefficients they 
also depend on the power spectrum index. The difference 
is that the 5^ parameters are averages of products of N p 
functions over the wave vectors, whereas the T, p param- 
eters are averages of the same products but over the or- 
thogonal part of the wave vectors. This latter problem is 
similar to the one that have been solved to get the 5^ . The 
only change is that the calculation has to be done in a 2D 
Fourier space instead of a 3D. To derive this series I follow 
the method developed by Bernardeau (1992, 1994b). 

2.3.1. The S p parameters without smoothing 

The first step of the calculation is to derive the generating 
function of the coefficient T, p when one makes the approxi- 
mation that W 2 (l) ?s 1. It will turn out that it corresponds 
to the case n = —2. For 3D calculations the correspond- 
ing case has been considered by Bernardeau (1992). The 
only difference with the 3D case is that the mean value 
of (ki.k 2 ) 2 1 (k\k 2 ) 2 is 1/2 instead of 1/3. The generat- 
ing function of the vertices (i.e., the monopole part of 



i\p[Ki , . . . , Kp)), ys\ i ), is men given uy une uinerenijiai 
equation, 



(27) 



with Qs(t) 



-t when r — ► 0, 



instead of equation (27) of Bernardeau (1992). I did not 
find any simple solution for this differential equation, but 
it can be shown that it corresponds to the equation de- 
scribing the "spherical" collapse in 2D, where Gs(t) is the 
density contrast and — r is the linear density contrast. It 
can also be shown that 



Gs(t) - 1 ~ r -(vi3-i)/2 w hen r- 
and that the form 



Gs(r) 



l+ T - 

v 



1 with v 



13- 1 



(28) 



(29) 



provides a good fit for the solution of equation (27). More 
precisely one can rigorously calculate the expansion of 
Gb(j) near r = and it reads 



12 



29 



79 



2085 



14 42 147 5096 
The generating function of the parameters T, p , 

i-y) p 



p P \ 



Si 



i, 



(30) 



(31) 



P =i 

is solution of the system, 

<p{y) = y + yGb[T{y)} + \ T iy) 2 
r{y) = -yO" s [T(y)]. 



(32) 



The system (32) is a standard result of tree summation 
(Bernardeau & Schaeffer 1992, Bernardeau 1992). Using 
the expansion (30) for the function Gs{j) one can calculate 
the expansion of ip(y) to get the first values of T, p . I found, 



E 3 (n 

S 4 (n = -2) 

S 5 (n = -2) 
S 6 (n = -2) 



2) = f; 



2540 ^ 
49 ' 
793; 
16370: 



(33) 



These results do not take into account the filtering ef- 
fects and are then independent of the shape of the power 
spectrum. They are actually exact for n = —2 since the 
integrales are then dominated by the k — ► limit where 
W 2 (xk6 )&l. 
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Fig. 1. The ratios E p as a function of the power law index, 7 
(solid lines) for p = 3, . . . , 6 and for an Einstein-de Sitter Uni- 
verse. The ratios S p (Eq. [1]) for 3D statistics are represented 
by dashed lines. 



2.3.2. 



The exact T, p series 



The derivation of the exact S p series has to take into 
account the geometrical dependences contained in the win- 
dow function (i.e. in the last section I assumed actually 
that !U 2 (|fci + . . . + fcp|) = W 2 (h) . ..W 2 (k p )). Calcula- 
tions made in 3D prove that the effects of the filtering are 
significative. Comparison of S3 in (33) and the result (22) 
proves also that it has to be taken into account. 

For a top-hat window function the properties (20) and 
(21) show that the smoothing corrections lead to analyti- 
cally simple corrections. It also demonstrates that we en- 
counter the same situation than for the derivation of the 
smoothing corrections in 3D with a top-hat window func- 
tion. Following the structure of the results obtained by 
Bernardeau (1994b) I expect that the smoothing correc- 
tions will be given by a Lagrangian-Eulerian space trans- 
formation. The generating function of the exact cumu- 
lants, ip(y), is then given by a system similar to (32) but 
with a generating function of vertices Gg (t) given by 



Gi(r) = Gs (r[l + Gi(r)] 



-(2+n)/4 



(34) 



±iie ueiivauion 01 une expressions 01 une nisi lew z^p is 
then quite straightforward and can be done easily with a 
mathematical symbolic package. 

I give here the expression of the first few parameters 



-'pi 



S, = — 



36 3(n + 2) 



7 2 
2540 



33 (n + 2) + 



21(n + 2) 2 



(35) 



49 v 4 

5 5 = 793 - 794 (n + 2) + 265 (n + 2) 2 

- 29.4 (n + 2) 3 ; 

5 6 = 16370 - 22511 (n + 2) + 11594 (n + 2) 2 

- 2650 (n + 2) 3 + 226.9 (n + 2) 4 ; 



In general there are no simple closed forms for the ratios 
S p . 

In Fig. 1, I present the dependence of the ratios S p as 
a function of the power law index. They are compared to 
the S p parameters (Eq. [1]) for the 3D statistics (dashed 
lines). It shows that they are significantly greater than 
their 3D counterparts. 

2.4. The £1 dependence of the s p parameters 

In the small angle limit the factorisation property (25) 
is always valid, whatever the value of £1. The parameters 
s p can then depend on SI through the parameters R p or 
the parameters S p . The SI dependence of R p appears only 
when the time evolution of the correlation functions is 
taken into account (that is for very deep surveys), since 
the radius-time relationship is SI dependent. But as in all 
the results that have been presented in this paper such 
a time dependence has been ignored, the SI dependence 
of R p will not be discussed furthermore and be assumed 
negligible. 

The SI dependence is then entirely contained in the S p 
parameters. To solve this problem in general we have to 
use the dynamics of the "2D spherical collapse" for any 
value of SI, which is not known. It is possible however to 
get, without too much difficulty, the value of S3 in the 
SI — ► limit. This is based on the result that in such a 
limit we have (see Bouchet et al. 1992), 



Po 



(SI = 0) 



exp[i(fci + k 2 )x 



d 3 fci d 3 fc 2 

(2tt) 3 / 2 (2tt) 3 / 2 

3 fci .k 2 1 
1 — - H — 

4 + k\ + 4 



Hki) 6(k 2 ) 



ki.k 2 

k 2 ki 



The resulting value of S3 is then 
=3(0 = 0) = ^-^. 



(36) 



(37) 



It demonstrates that for S3 the SI dependence is extremely 
weak: between SI = 1 and SI — ► the variation of S3 is 



ouiy 01 une oruer 01 0/0. iviore precisely loiiowing rjouciieu 
et al. (1992) one can use an analytical fit to describe the 
SI dependence of S3 valid for 0.05 £ SI ^ 3., 



S 3 (S1) = f + ^ (J2-2/63 _ ! 



3(n + 2) 



(38) 



A similar weak SI dependence is expected for the higher 
order parameters S p (as it is the case for the S p coeffi- 
cients). 

3. Validity of the small angle approximation 

The expressions (22-23) of the angular skewness and the 
higher order coefficients (26, 35) rely on the approxima- 
tion (15) which corresponds to the small angle limit. The 
validity domain of this approximation is however ques- 
tionable. This investigation is all the more crucial that in 
principle the quasilinear regimes applies at large physical 
scales which would correspond to rather large smoothing 
angles. 

To test the validity of this approximation I computed 
the S3 parameter from direct Monte-Carlo numerical inte- 
gration for different smoothing angles and sample charac- 
teristics. 

3.1. Model for galaxy clusters 

To start with I used a very simple model given by, 



^ m0d » = ( - 



(39) 



with 7 = 1.5, 1.7 or 1.9, and ro = 20/i -1 Mpc. The selec- 
tion function is simply assumed to be a Heavyside func- 
tion, 



F model(1 \r) = 1 for r < D 



F 



model(l) 



(r) = for r > D 



(40a) 



(406) 



with D r; 600/i -1 Mpc. This could actually be considered 
as a good parametrization of the properties of galaxy clus- 
ter samples (Toth et al. 1989). The other advantage of 
such a choice is that it allows to test the Monte-Carlo inte- 
gration against analytical results in two cases: in the small 
angle limit, and for #0 = 180°, the latter corresponding to 
the filtering with a spherical top-hat window function, (it 
has obviously no practical interest for angular surveys!) 

The results for the mean 2-point angular function are 
presented in Fig. 2. The solid lines show the results of 
Monte-Carlo integrations for different angles and different 
values of 7. The dashed lines correspond to the small angle 
Limber approximations, 
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Fig. 2. The mean 2-point angular correlation function as a 
function of the smoothing angle for the model (1) (39, 40). 
The thick solid line is for 7 = 1.7 and the thin solid lines for 
7 = 1.5 and 7 = 1.9. The dashed lines correspond to the small 
angle limit (41) and the circle to the #0 = 180° case (eq. [42]). 
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Fig. 3. The value of S3 as a function of the smoothing angle for 
the model (1) (39,40). The notations are the same than in Fig. 2 
where the dashed lines correspond here to the small angle limit 
value, (44) and the circles to the expression, S3 = 34/7 — 7. 

and the circle to the exact integration corresponding to 
the case of the spherical top-hat window function 



^(180°)= 7 



9 2 3 "T 



( 6 4) ( 7 - 3) 



(42) 



The mean three-point angular correlation function has 
been obtained using the expression of the three-point cor- 
relation function in real space given by the second order 
perturbation theory (Fry 1984), 



1 



cos(0i) — + — 
V r 3i r 2 i 



3 - 27 + 7 2 cos 2 I 



(43) 



7 (3- 7 ) 2 
+ cyc.(l,2,3) 

where 9\ is the angle between the vectors r 2 — r\ and 
rs — r\ and r 2 \ and rz\ are the lengths of these vectors. 
For the parameters of this model the equations (22, 23) 
give, for an Einstein-de Sitter Universe 1 , 



model(l) 



(Oo - 0) 



(6-T) 2 



93 3 
14 ~ 2 7 



(44) 



3(9 - 2 7 ) 

The resulting values for S3 (= ZZT 3 /ZZT| ) are presented in 
Fig. 3 as a function of #0 for the different values of 7. The 
notations are the same than in Fig. 2. It has to be noticed 
that the small angle limit is reached only for very small 
angle, e.g. significantly smaller than 1°. The transition 
between the small angle result and the #0 = 180° case 
happens for 6 ps 1°. 

It implies that the analytical results of the previous 
part cannot be applied exactly to the observations in the 
1° ~ 5° range for such a model. However this result sug- 
gests that the exact values are intermediate between the 
small angle limits and the value of the S p parameters (Fig. 
1). It thus gives a plausible range where the exact values 
of the s p parameters are expected to be. 

3.2. Model for galaxies 
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Fig. 4. The mean 2-point angular correlation function as a 
function of the smoothing angle for the model (2), eq. (45), of 
the selection function and for 7 = 1.7. The triangles correspond 
to the observations in the APM galaxy survey (Gaztanaga 
1994) 



1 In the limit SI 
(27/4- 3 7 /2). 



the factor (93/14 - 37/2) becomes 
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Fig. 5. The value of S3 as a function of the smoothing angle 
for the model (2). The dashed lines correspond to the small 
angle limits; the thin lines to the Q — > case; the thick lines to 
the Q = 1 case. The triangles correspond to the observations 
in the APM galaxy survey (Gaztanaga 1994) 
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the expected validity domain of the perturbation theory 
is for #0 ^ 1° (which corresponds to a real scale greater 
5/i _1 Mpc in which perturbation theory has been proved 
valid, Bernardeau 1994b, Juszkiewicz et al. 1994, Baugh 
et al. 1994). 

When they are compared to the observations the per- 
turbation theory predictions for S3 in Fig. 5 are found to 
be greater (4.6 instead of 3.8 ± 0.2 for 6 « 2°) in the 
1 — 5° range. This excess is only about 20% but leads to 
a significative disagreement between perturbation theory 
and the observations, according to the errorbars claimed 
by Gaztanaga, and contrary to the conclusions of the anal- 
ysis made by Frieman & Gaztanaga (1994). However, it is 
important to have in mind that the perturbation theory 
predictions obtained in the previous section neglected ef- 
fects such as the redshift evolution of the correlation func- 
tions, the departure of the two-point correlation function 
from a power law bahaviour, errorbars for the determina- 
tion of the 2-point correlation function. These effects may 
change the results. Thus more precise studies have to be 
undertaken before any definitive conclusion can be given. 



The second model I considered has been built to mimic 
the properties of the galaxy angular surveys. The selection 
function I used is then of the form 



Table 1. The s p parameters: theory and observation 



^model(2)( r ) ^ r -0. 



exp 



D 



(45) 



For such a model the small scale limit of the coefficient S3 
is, for an Einstein-de Sitter Universe, 



n model(2) , 



32 r(5/4) 
' 27 3 3 / 4 



93 



r(f-7) 
r2 (V) 



14 2 7 



(46) 



For the 2-point correlation function I assumed a power law 
behaviour with 7 = 1.7 and ro = 5/i _1 Mpc. It reproduces 
reasonably well the observed angular mean two-point cor- 
relation function (Fig. 4.) in the APM galaxy survey for 
#0 £ 2°. The resulting S3 parameters are given in Fig. 5. It 
can be seen that the small angle approximation is slightly 
better than for the previous case. It is not too surprising 
since the boundaries are in this case smoother, so that 
the small angle approximation (which basically neglects 
the effects of boundaries for 2 points out of 3) is expected 
to be more accurate. In such a case the values found at 
the 1° scale are in reasonable agreement with the small 
scale limit, although slightly lower. 



4. Comparison with observations and conclusions 

The measures of such s p coefficients have been made in 
the APM survey by Gaztanaga (1994) up to p = 9 in the 



p 


R p 


sp(e^o) 


obs: s p (9 > 1°) 


3 


1.19 


4.87 


3.81 ±0.07 


4 


1.52 


47.6 


32.5 ±4.2 


5 


2.00 


713.5 


384 ±62 


6 


2.71 


14730 


3260 ± 1340 



In table 1 I give the expected values of the s p (p = 3, 6) 
coefficients in the small angle limit using the parameters 
Rp of Gaztanaga (1994) for 7 = 1.7. The actual values 
predicted by the perturbation theory are expected to be 
close to these numbers when #0 ~ 1°. The observations 
are seen to be systematically smaller than the predictions, 
similarly to the p = 3 case. These predictions, however, 
rely entirely on the small angle approximation. 

If we admit that the discrepancy between perturba- 
tion theory results and observations is real, what would 
it mean? One could argue that perturbation theory does 
not apply, that nonlinear corrections affect the result. I 
think that this is unlikely to be the case since the main 
contribution to the averaged angular correlation functions 
comes from galaxies being at a distance of about 60 D 
which for the 1° scale is large enough for the perturba- 
tion theory to apply. Note also that the nonlinear correc- 
tions are rather expected to increase the values of the s p 
parameters and thus amplify the discrepancy. A more seri- 
ous concern would be the hypothesis that the galaxies are 
not biased. In this case it does not mean only that the so 
called b parameter should be 1, more constrainly it means 
that the density contrast observed in Galaxy counts re- 



piouuces une iiu nmieui aspecus 01 une maimer uensiuy con- 
trast (see for instance the discussion by Fry & Gaztanaga 
1993). Even if most of the mass of the universe is in galax- 
ies (possibly in haloes), the simple fact that there might 
be a segregation in luminosity for the correlation func- 
tion (Hamilton 1988, Park et al. 1994, Loveday et al. 1994 
but see also Alimi, Valls-Gabaud & Blanchard 1988 for 
a different conclusion) is likely to affect the whole statis- 
tics. The basic reason is that the various contributions are 
not properly weighted (i.e. by mass) when the correlation 
functions are measured. Bernardeau & Schaeffer (1992) 
explored the consequences of such effects on the two- and 
three-point correlation functions in the frame of a gravita- 
tionally induced luminosity segregation. They found sig- 
nificative changes for the value of S3. In particular, for the 
limit of rare objects (either very massive clusters or very 
bright galaxies) the value of S3 is expected to be close to 3 
for any cosmological models. This limit is not likely to be 
relevant for galaxy surveys but may provide a good model 
for the observed galaxy cluster correlation functions. 

In any case, even if the agreement between the obser- 
vations and the perturbation theory results is not perfect, 
it supports the idea that the present large-scale struc- 
tures have been built out by the groth of initial Gaussian 
fluctuations in a gravitational instability scheme. The s p 
coefficients provide then a unique and powerful way to 
test the large-scale structure formation models for gravi- 
tational instability scenarios. 
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